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I. INTRODUCTION 



The goal of this thesis is to study the combustion 
behavior of porous graphite . The mathematical model devel- 
oped by VATIKIOTIS (1) is a transient, one-dimensional model 
of a porous medium consisting of spherical particles , or 
cylindrical fibers of graphite . 

The model covers several aspects : mass transfer, heat 
transfer, combustion , airflow , and temperature dependency 
of thermophysical properties. Geometric parameters of the 
porous medium include porosity (p) and medium thickness (L). 
The porous medium has interconnected pores which permit air 
flow through the medium. This air flow was modeled by the 
continuity equation and Darcy’s Law . The porous medium is 
characterized by a number of parameters including porosity , 
specific internal area (Z), tortuosity (r) and permeability 
(m) . These parameters are defined as follows. Given a unit 
square cell of dimension D, and particle diameter d, we have 



= 1 - 7T / 4 (d/D) 2 



( 1 . 1 ) 



z = 1/2 v d 2 /D 3 



( 1 . 2 ) 



= 1.4 



(1.3) 



m = Sp 3 /z 2 



(1.4) 



Porosity,* p, is defined as the ratio of void volume per 
volume of unit all . The specific internal area, Z, is the 
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ratio of internal surface area to bulk volume . The 
tortuosity is defined as the ratio of the length of the 
flowpath of a fluid particle to the straight line distance. 

The particle's size decreases as the carbon is consumed 
and all geometric properties which depend on particle diam- 
eter are functions of time and position . This model 
considers that the carbon matrix remains rigid as the 
particle diameter decreases and thus porosity increases with 
combustion . 

The magnitude of REYNOLDS number defines the motion of 
the fluid which can be molecular , viscous or inertial. In 
the case of porous media , the flow is dominated by viscous 
and inertial effects. Darcy's Law for the range of REYNOLD'S 
number where viscous effects dominate is the equation that 
governs the fluid flow in porous media. Two velocities which 
distinguish the flow through the medium are filter velocity 
Q, and pore velocity u. Filter velocity is given by : 

Q = - mfn ( dP/dx ) (1.5) 

where u is the viscosity of the fluid and P is the pressure. 

The pore velocity u is given by : 

U = Q/p (1.6) 

The hypothesis of the Dupuit-Forcheimer assumption is 
that the local pore velocity is greater than the filter 
velocity. The actual velocity in a pore is a function of 
position within the pore . The Dupuit-Forcheimer relation 
defines an average velocity in a pore. 

The model of combustion adopted was the model of N. N. 
Semenov (2). The relation of reaction rate to temperature 
and oxygen concentration , and the interaction of heat gener- 
ation and heat transfer are fundamental to this model . If 
the heat transfer dominates, extinction will occur. If the 
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heat generation dominates then combustion will result. The 
reaction rate equation was adopted from Arrhenius's Law for 
a simple reaction where the rate depends on the concentra- 
tion of the reactants and not on the products. The heat 
generated by an exothermal reaction is obtained by multi- 
plying the reaction rate by the heat of combustion . 

A study of the curve of heat generation versus tempera- 
ture presented by Frank-Kamenetskii distinguishes two phases 
during combustion. During the initial kinetic phase , the 
rate of reaction and the temperature are lowest. During This 
phase , there is an excess supply of oxygen. In the kinetic 
regime the reation rate increases exponentially with 
increasing temperature. The kinetic regime is followed by 
the diffusion regime. This phase is distinguished by high 
temperatures and reaction rates , and the reaction is 
limited by a lack of oxygen. 

It is assumed that the chemical reaction for the combus- 
tion process produces carbon monoxide and carbon dioxide. 
The ratio of the mass rates of carbon monoxide to carbon 
dioxide depends on the temperature. Increasing temperature 
results in an increase of this CO to CO 2 ratio. 

The model of heat transfer includes three basic mecha- 
nisms : convection , conduction and radiation . The heat 
transfer equation used by VATIKIOTIS (1) also includes a 
heat generation term due to combustion. All thermophysical 
properties were treated as temperatures dependent. These 
properties include viscosity, conductivity and density. 
Energy balance equations were constructed for control 
volumes of the porous medium and the air flowing through the 
porous medium . 

A third field equation of the model was obtained by a 
mass balance of oxygen . This equation includes molecular 
diffusion, convective transport and a sink term to account 
for the depletion of oxygen due to combustion. 
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These equations resulted in a system of coupled, nonli- 
near transient field equation . These field equations, 
together with boundary and initial conditions , define the 
problem. They were solved numerically by using the Galerkin 
finite element method. Details of the formulation are 
presented in reference 1 . 
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II. EFFECT OF HEAT INPUT ON COMBUSTION 



A. INTRODUCTION AND DESCRIPTION 

In this section the effect of heat input on the combus- 
tion behavior of porous medium was investigated. The porous 
medium is subjected to a heat source on its surface . Two 
parameters of main importance are considered: The magnitude 
of the heat flux of the source to the porous medium (SQ), 
and the time that the heat flux remains active (TQ) . 

To do this analysis, all parameters of the system are 
fixed except the time period of heating (TQ) and the amount 
of heat flux (SQ) . 

The purpose of this section is to determine the relation 
between SQ and TQ which will result in combustion. This is 
achieved as follows: 

For a specific value of heat flux SQ, the program is run 
with different values of TQ until combustion occurs at (TQ) C 
and extinction occurs at (TQ) e , where5 is a small time 
period . Here we took<5 seconds . 

TQ C = TQ e + 5 (2.1) 
Four cases of heat flux were studied in this section: 



CASE 


II. 1 


SQ = 


18000 


BTU/ft . sq 


hr 


CASE 


II. 2 


SQ = 


20000 


BTU/ft . sq 


hr 


CASE 


II. 3 


SQ = 


30000 


BTU/ft.sq 


hr 


CASE 


II . 4 


SQ = 


40000 


BTU/ft . sq 


hr 



The fixed parameters used in all cases are: 

Ambient temperature =80 F 
Ambient temperature = 2117 Ib/ft.sq 
Tortuosty = 1.400 
Filament diameter = 0.0004167 ft 
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Thickness of matrix laminate = 0.0004167 ft 
Thickness of porous medium = 0.02083 ft 
Gas constant for air = 53.34 lbf ft/lbm R 
Conductivity of filament = 86 BTU/lbm H.F. 

Specific heat of filament = 0.703 lbm/Cf 
Emissivity of filament = 0.90 

Shape factor for int.HF.XFER coefficient = 1.00 
Characteristc length of medium = 1 inch 
Heat of reaction = 14090 BTU/Cf 
Reaction order = 0.50 

Stochiometric ratio (fuel/air) = 0.375 

Reaction coefficient = 2065000 lbm/Cf H 

Activation energy coefficient = 28840 deg R 

Pressure differential across thickness = 50 lb/ft. sq 

Initial carbon and air temperature =80 F 

Initial uniform oxygen concentrat ion = 0.0172 lbm/ft.cu 

The equations used in this program permit three sets of 
boundary conditions that are approximations to physical 
situations. The boundary conditions selected are the best 
approximations of a one dimensional model. The set of 
boundary conditions used in this case are typical of thermal 
flow reactor without radiation from the boundary surfaces . 
for the problems of this investigation, the following 
boundary conditions were used : 

- No heat transfer from either end of the porous medium, 
dt/dx =0 at X/L = 0 and X/L = 1 ( 2.2 ) 



- Danckwerts' boundary conditions : Conditions for the 
air temperature and oxygen concentration were assumed. This 
means : 
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K a 3T a /3x = X/L = 0 (2.3) 

3T a /dx =0 X/L = 1 (2.4) 

deft / c)x - u ~ *P * os') X/L - 0 (2.5) 

(K e +K r ) aT c /ax = -a e(T c 4 - T 4 ^ X/L = 1 (2.6) 

K a a T a /ax = p. a c a u(T a - Too) X/L = 1 (2.7) 

aT a /ax = 0 X/L = 1 (2.8) 

D e a<£/ax = u (0 - ^ X/L = 1 (2.9) 

B. PROCEDURE 

For each case , the amount of heat flux (SQ) was fixed 
and several values of shut off time (TQ) were chosen until 
values of TQ C and TQ e were found such that : 

TQ C * TQ e (2.10) 

TQ C = TQ e + § (2.11) 

where 5 is a small number. 
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C. RESULTS 



The results for each case were obtained in both numer- 
ical and graphical form . The behavior of carbon tempera- 
ture, oxygen concentration and reaction rate for combustion 
and extinction are shown as surfaces versus time and posi- 
tion The program also calculates the air temperature in the 
porous medium. However the air temperature surfaces are not 
shown as they are quite similar , in most cases to carbon 
temperature . 

There are three surface graphics for each heat flux 
case : 

-Carbon temperature versus position and time 

-Oxygen concentration versus position and time 

-Reaction rate versus position and time 

In adition , there are two curves for each heat flux 
(SQ) , one showing the variation of carbon temperature 
versus time for combustion and the other showing the varia- 
tion of carbon temperature for extinction . These is two 
curves are for the point x/L=0. (i.e. where the air enters 

the porous medium) . 

1. CASE II-l SQ = 18000 BTU/ft.sq hr 

The values for shut off time bounding extinction and 
combustion are : 



for 


combustion 


TQ C 


for 


extinction 


TQ e 



46 seconds 

47 seconds 



During the transient 
varied . for each analysis , 
of some of these parameters 
character of the particular 
their values are given in Tab 
The graphical results 



analyses , system parameters 
the maximum and minimum values 
provide some insight into the 
case . These parameters and 
Le I 

are shown in Figures 2.1 to 



2.6 
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u- 



? pnr 




Figure 2.1 Temperature vs X/L and 
for SQ = 18000 BTU/ft.sq hr and TQ = 46 

Tfr^rpsTijpr 




Figure 2.2 Temperature vs X/L and 
for SQ = 18000 BTU/ft.sq hr and TQ = 45 



t ime 
seconds 



time 

seconds 
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OXYGEN GONE. SUREFCE FROM GRRF3E 




Figure 2.3 Oxygen concentration vs X/L and time 
for SQ = 18000 BTU/ft.sq hr and TQ = 46 seconds. 

on GEN CONG. 5 'J? rc CE FROM u 3 Fr it 




Figure 2.4 Oxygen concentration vs X/L and time 
for SQ = 18000 BTU/ft.sq hr and TQ = 45 seconds. 



PlRCTION PR?: SURPFCC r Pj w 




Figure 2.5 Reaction rate vs X/L and time 
for SQ =18000 BTU/ft.sq hr and TQ = 46 seconds. 

peirction prt: 5u- r R:: r:? c~ :== r 3: 




Figure 2.6 Reaction rate vs X/L and time 
for SQ = 18000 BTU/ft.sq hr and TQ = 45 seconds. 



25 



GRAPHITE TEMPERATURE (K) 

1000 1X00 1500 1700 SOOO 8X0 0 1600 8700 3000 



TIME VERSUS GRAPHITE TEMPERATURE 

o 




Figure 2.7 Combustion and extinction carbon temperature 

for position X/L = 0.. 
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2. CASE II-2 SQ = 20000 BTU/ft.sg hr 

For this case the heat flux duration times for the 
conditions of extinction and combustion are : 

for combustion TQ C = 35 seconds 

for extinction TQ e = 34 seconds 

During the transient history of this case several 
parameters varied. The maximum and minimum of these changing 
variables are in the Table II . It is necessary to emphasize 
that these extrema results do not occur at the same time or 
same position (X/L) . They are maximum and minimum values 
obtained at different times and position during the tran- 
sient process. These values show the characteristc behavior 
of this case . 

The behavior of carbon temperature , oxygen concen- 
tration and reaction rate are shown in Figures 2.8 to 2.13 . 
These surfaces are an important characteristic of this case 
because they show the evolution of the combustion process or 
extinction process for each position during the transient 
t ime . 

3 • CASE II- 3 SQ = 30000 BTU/ft.sq hr 

The values of heating time that result in combustion 
and extinction for this value of heat flux are : 



for 


combustion 


TQ C 


= 15 seconds 


for 


extinction 


TQ e 


= 14 seconds 



During the transient analysis , some variation of 
values occured to several parameters. The limits of varia- 
tion , maximum and minimum values , are shown in Table III . 
The process of combustion and extinction have particular 
characteristics that determine the variation of the values 
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TOIPERPTUPC SURFRCC fRC“ G D R r ZE 




Figure 2.3 Temperature vs X/L and 
for SQ = 20000 BTU/ft.sq hr and TQ = 35 

TEfiPCRFlTURL SURFRCl FROM GRRF3C 




Figure 2.9 Temperature vs X/L and 
for SQ = 20000 BTU/ft.sq hr and TQ = 34 
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OXYGEN CONC. SUPEECE ERG" GREE3E 




Figure 2.10 Oxygen concentration vs X/L and time 
for SQ = 20000 BTU/ft.sq hr and TQ = 35 seconds. 

OXYGEN CONC. SURFRCE r R0^ GREECE 




Figure 2.11 Oxygen concentration vs X/L and time 
for SQ = 20000 BTU/ft.sq hr and TQ = 34 seconds. 
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REACTION RATE SURFACE r RC'M CRA- - 3F 




Figure 2.12 Reaction rate vs X/L and time 
for SQ = 20000 BTU/ft.sq hr and TQ = 35 seconds. 

REACTION RATE SURFACE FROM GRA r 3E 




Figure 2.13 Reaction rate vs X/L and time 
for SQ =20000 BTU/ft.sq hr and TQ = 34 seconds. 
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TIME VERSUS GRAPHITE TEMPERATURE 

O 




Figure 2.14 Extinction and combustion carbon temperature 

for position X/L = 0. . 
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of these parameters. The maximum and minimum values give a 
general idea of system behavior . These limiting values do 
not occur at same position (X/L) or same time . 

A better idea of the development of temperature , 
oxygen concentration and reaction rate can be obtained from 
Figures 2.15 to 2.20 These surfaces show the values of these 
parameters for each point (X/L) during the transient time . 

Figure 2.14 shows that at X/L = 0. the transient 
behavior for extinction and combustion are the same until 
TQ e = 14 seconds , the time the heat flux is terminated . 
After this time the temperature for the extinction case 
decreases to ambient temperature ; while for the combustion 
case , the temperature increases reaching a peak temperature 
of 3200 F at 28 seconds and thereafter decreases to an equi- 
librium combustion temperature of about 2500 F after about 
50 seconds (see Figure 2.14). 



4. CASE II-4 SQ = 40000 BTU/ft.sq hr 

The heat flux of this case produced combustion and 
extinction of the porous medium at the following times : 

for combustion TQ C = 9 seconds 

for extinction TQ e = 8 seconds 

During the transient process , the output of the 
program showed parameters th at didn't change with time and 
others that did . The variation of these parameter's values 
show the development of the combustion and extinction 
phenomenon for this value of SQ . The maximum and minimum 
values of several parameters are shown in Table IV . Again 
these values do not correspond to the same time . 

The characteristics of this case , carbon tempera- 
ture , oxygen concentration and reaction rate , as function 
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TEMPERATURE SURFACE FROM GRAF3E 




Figure 2.15 Temperature vs X/L and time 



for SQ = 30000 BTU/ft.sq hr and TQ = 15 seconds 



TEMPERATURE SURFACE FROM GRAF3E 




Figure 2.16 Temperature vs X/L and time 
for SQ = 30000 BTU/ft.sq hr and TQ = 14 seconds 
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OXYGEN CONG. SURFACE FROM GRAF3E 




Figure 2.17 Oxigen concentration vs X/L and time 
for SQ =30000 BTU/ft.sq hr and TQ = 15 seconds. 

OXYGEN CONC. SURFACE FRO* GRAF3E 




Figure 2.18 Oxigen concentration vs X/L and time 
for SQ = 30000 BTU/ft.sq hr , and TQ = 14 seconds. 
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RERCT10N RATE SURFACE FROM GRAF3E 



Figure 
for SQ = 



Figure 
for SQ = 




2.19 Reaction rate vs X/L and time 
30000 BTU/ft.sq hr and TQ = 15 seconds 



REACTION RATE SURFACE FRO" C-RFF3E 




2.20 Reaction rate vs X/L and time 
30000 BTU/ft.sq hr and TQ = 14 seconds 
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GRAPHITE TEMPERATURE (K) 

*60 600 760 10O0 1*60 1600 1760 *000 **60 *600 *760 *000 3260 



TIME VERSUS GRAPHITE TEMPERATURE 




2.21 Extinction and combustion carbon temperature 
for position X/L = 0.. 
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of time and position (X/L) are shown as graphic surfaces in 
Figures 2.22 to 2.27 . These graphical results give an idea 
about the behavior of each point of the transient history . 

For X/L = 0 , both the combustion and extinction 

cases show similar behavior to this time equal to about 8 
seconds , when the heat flux is terminated . After this 
time, for the extinction case , the temperature decreases to 
ambient temperature . For the combustion case , the temper- 
ature at X/L = 0 increased to about 3300 F , and thereafter 
decreased to around 2500 F (temperature of combustion) 
Both the extinction and combustion cases achieve new equi- 
librium temperatures , about 2500 F for the combustion cases 
and ambient temperature for the extinction cases , after 35 
seconds (see Figures 2.28 ). 

D . SUMMARY 

Here , some observations are made about the effects of 
heat flux input and duration of heat flux on the combustion 
process . 

1 . Power Relation 

For each case , a pair of values (SQ,(TQ) C ) was 
obtained . A plot of these point is shown in Figures 2.29 
and 2.30 . It is observed that the relation between SQ and 

TQ of these graphics yields the approximate power relation 
( equation 2.12) . 



2 . Relation of Temperature and Oxygen Concentration 

During the initial heating of the porous medium , 
the temperature increases while the oxygen is being 
defleted. After the heat flux is removed , either extinction 
or combustion will proceed. If extinction occurs , the 
carbon temperature decreases to ambient temperature while 
the oxygen is being restored to ambient concentreition . 

If combustion occurs , then the oxygen inside the 
medium goes to zero and remains there . The only oxygen 
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TEMPERATURE SURFACE FROM 3RRF3E 




Figure 2.22 Temperature vs X/L and 
for SQ = 40000 BTU/ft.sq hr and TQ = 9 

TEM°EPRTURE 5U 5r ACE FRC M 3 C "3E 




Figure 2.23 Temperature vs X/L and 
for SQ = 40000 BTU/ft.sq hr and TQ = 8 



t ime 
econds . 



time 

seconds . 
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OXYGEN CONC. SURFRCE FROM GRRF3E 



Figure 2 
for SQ 



Figure 2 
for SQ 




= 40000 BTU/ft.sq hr and TQ = 9 seconds. 
OXYGEN CONC. SURFRCE "ROM GR C "3E 




.25 Oxigen concentration vs X/L and time 
= 40000 BTU/ft.sq hr and TQ = 8 seconds. 
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REACTION RATE SURFACE FROM 3RAF3E 




F igure 
for SQ = 



2.26 Reaction rate vs X/L and time 
40000 BTU/ft.sq hr and TQ = 9 seconds. 



REACTION RATE SURFACE FROM GRBF3E 




Figure 2.27 Reaction rate vs X/L and time 
for SQ = 40000 BTU/ft.sq hr and TQ = 8 seconds 
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GRAPHITE TEMPERATURE (K) 
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TIME VERSUS GRAPHITE TEMPERATURE 




Figure 2.28 Extinction and combustion carbon temperature 

for position X/L = 1. 
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3.4 10 10 / ( SQ) 2 . 09 



( 2 . 12 ) 



TQ C = 

present is at the entrance at X/L = 0 . The temperature 
continues to increase to ' a maximum , and thereafter the 
temperature decreases . Eventually the entire medium would 
achieve a uniform combustion temperature . 

3 . Equilibrium Temperature for Combustion 

In all combustion cases , the behavior was similar . 
During the initial combustion period , the temperature of 
the medium at X/L = 0 rose between 2000 F (CASE II-l) to 

3400F (CASE II-4) . At this time , on the other side of the 
medium at X/L =1. , the temperature of the medium was 

approximately 1500 F . This during thus early transient 
stage there is aAT of 1300 F (for CASE II-l) and 1900 F 
(for CASE II-4) across the medium . As time proceeds the 
temperature at X/L = 0 decreased , while at X/L =1. , the 

temperature increased . That is , with progressing time , 
the medium moves towards a uniform temperature . In the 
present case , this equilibrium temperature of combustion is 
around 2400 F to 2500 F for all cases . 

The mechanism of equilibrium of the temperature 
during combustion is due to the heat transfer by air . As 
the air flows through the porous medium it transports the 
heat from the hot front part of the medium to the cooler 
back part , until all points of the medium have the same 
temperature . 
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m i shut or r time; hcptcrj 



SQ VERSUS Td 




MINIMUM TO TOR COMBUSTION 



Figure 2.29 Heat flux versus combustion temperature (TQ C ) 

rectangular plot. 
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TQ (SHUT OFF TIME HEATER) 



SQ VERSUS TQ 



o 




Figure 2.30 Heat flux (SQ) vs combustion temperature (TQ C ) 

Log log plot. 
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III. EFFECT OF THICKNESS ON COMBUSTION 



A. INTRODUCTION 



The porous medium has several geometric parameters that 
can have direct influence on the combustion process. The 
effect of the thickness parameter on system behavior is 
analysed in this section. All other parameters are fixed 
except the thickness of the medium. 

The purpose of this section is to determine the relation 
of the thickness on the minimum temperature to start combus- 
tion (i.e., the ignition temperature). 

For a specified value of thickness, the program is run 
with different values of initial temperature. This initial 
temperature is equal at all points of the medium. The lowest 
value of initial temperature that results in combustion is 
the ignition temperature Tc. Extinction occurs at Tc-<5 , 
whereS is a small increment of temperature. The value of<5 
used in this parametric analysis was 10 F . 

Seven cases with different thicknesses were studied in 
this section: 



CASE III-l 
CASE III- 2 
CASE III- 3 
CASE III-4 
CASE III- 5 
CASE III- 6 
CASE III-7 



THICKNESS 

THICKNESS 

THICKNESS 

THICKNESS 

THICKNESS 

THICKNESS 

THICKNESS 



0,25 inches 
0,50 inches 
0,75 inches 

1.00 inches 

2.00 inches 

4.00 inches 

6.00 inches 



Each case is characterised by its thickness. The fixed 



parameters used in all cases are: 



- Ambient temperature = 80 F 

- Ambient pressure = 2117 lb/ft. sq. 

- Tortuosity = 1.400 

- Filament diameter = 0.0004167 ft. 
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- Thickness of matrix laminate = 0.0004167 ft. 

- Thickness of porous medium = 0.02083 ft. 

- Gas constant for air = 53.34 lbf.ft/lbm.R 

- Conductivity of filament = 86 BTU/lbf HF 

- Specific heat of filament = 0.703 lbm/cf. 

- Emissivity of filament = 0.90 

- Shape factor for int.HF.XFER coefficient = 1.00 

- Heat of reaction = 14090 BTU/CF 

- Reaction order = 0.50 

- Stoichiometric ratio (full/air) = 0.375 

- Reaction rate coefficient = 2065000 lbm/CF hr 

- Activation energy coefficient = 28.840 deg R 

- Pressure differential across thickness = 50 lb/ft 

- Initial uniform oxygen concent rat ion = 0.0172 lbm/ft 

B . PROCEDURE 

For each case, the characteristic length of the medium 
(thickness) was fixed and values of the uniform initial 
temperature were chosen until the value of temperature for 
combustion ( T c ) and temperature for extinction (T e ) were 
approximately equal (T c « T e ) or (T c = T e + 10) (F) . 

C. RESULTS 

For each case, the results were obtained in both numer- 
ical and graphical form. The behavior of carbon temperature, 
oxygen concent rat ion and reaction rate for combustion and 
extinction are shown as surfaces versus time and position. 

The numerical results from the program show the varia- 
tion of several parameter values during the transient anal- 
ysis. The air temperature inside the porous medium, the 
griphite temperature and the oxygen concentration are 
presented at each position (X/L) at intervals of time. Other 
system properties are also given . These properties charac- 
terize the process of combustion and extinction for 
different thicknesses. The maximum and minimum do not occur 
at the same point or at the same time. 
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The graphical results are given in the form of surfaces. 
These surfaces show the behavior of carbon temperature, 
oxygen concentration and reaction rate versus time and posi- 
tion (X/L) for the combustion and extinction processes. 

1. CASE III-l Thickness = 0.25 inches 

The values of the initial uniform carbon temperature 
bounding the extinction and combustion processes are: 



for combustion Tc = 1100 F 
for extinction Te = 1090 F 



The development with time of the behavior of the 
graphite temperature, oxygen concentration and reaction rate 
during the combustion and extinction processes are shown in 
Figures 3.1 to 3.6 The transient parametric analysis provide 
some insight into the characterers of this particular case. 
The maximum and minimum values of these parameters are given 
in Table V . 

2 . CASE III-2 Thickness = 0.50 inches 

The combustion and extinction process, in this case, 
start with the following uniform carbon temperatures: 



for combustion Tc = 960 F 
for extinction Te = 950 F 



For this case, the graphical results for graphite 
temperature, oxygen concentration and reaction rate are 
shown in Figures 3.7 to 3.12 . The characteristics of the 
changing parameters in this case, are shown by the maximum 
and minimum values in Table VI . 
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TEMPERATURE SURFACE FROM GRAF3E 




Figure 3 . 1 Tsmperature vs X/L and time 
for thickness = 0.25 in 
Initial carbon temperature = 1100 F. 

TEMPERATURE SURFACE FROM GRAF3E 




Figure 3.2 Temperature vs X/L and time 
for thickness = 0.25 in 
Initial carbon temperature = 1090 F. 
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OXYGEN CONC. SURFRCE FROn GRRF3E 



Figure 



Figure 




3.3 Oxygen concentration vs X/L and time 
for thickness = 0.25 in 
Initial carbon temperature = 1100 F. 

OXYGEN CONC. SURFRCE FROM GRRF3E ; 




3.4 Oxygen concentration vs X/L arid time 
for thickness = 0.25 in 
Initial carbon temperature = 1090 F. 
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REACTION RATE SURFACE FROM GRAF3E 




Figure" 3.5 Reaction rate vs X/L and time 
for thickness = 0.25 in 
Initial carbon temperature = 1100 F. 

REACTION RATE SURFACE FROM GRAF3L 




Figure 3.6 Reaction rate vs X/L and time 
for thickness = 0.25 in 
Initial carbon temperature = 1090 F. 
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TEMPERATURE SURFACE FROM GRAF3E 




Figure 3.7 Temperature vs X/L and time 
for thickness = 0.50 in 
Initial carbon temperature = 960 F. 

TEMPERATURE SUR r PCE FROM GRAF3E 




Figure 3.8 Temperature vs X/L and time 
for thickness = 0.50 in 
Initial carbon temperature = 950 F. 
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OXTGEN CONC. SURFACE FROM GRRF3E 




F igure 



3.9 Oxygen concentration 
for thickness = 0.50 
Initial carbon temperature 



vs X/L and time 
in 

= 960 F. 



GXTGEh CONC. SURF-CE 



-“ROM GRRF3E 




Figure 3.10 , Oxygen concentration vs X/L and time 
for thickness = 0.50 in _ 

Initial carbon temperature = 950 F. 
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Ruction rate surface: 



p ROf1 GRAF3C 




figure 3.11 



Ini 



Reaction rate vs 
for thickness = 0.50 
tial carbon temperature 



R/L and time 
in 

= 960 F. 



Reaction rate surface: 



ERO M GRAF3C 



Figu 




re 3.12 



Reaction rate vs X/L and time 
for thickness = 0.50 in 
Initial carbon temperature = 950 F. 



38 



TABLE VI 
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3. CASE III - 3 Thickness 



0.75 inches 



During this case the minimum uniform temperature 
that starts the combustion process and the maximum uniform 
temperature for extinction are: 

for combustion Tc = 900 F 
for extinction Te = 890 F 

During the transient analysis, the system parameters 
varied. For this case , the maximum and minimum values of 
some of these parameters provide some insight into the char- 
acter of the case. These parameters and their extreme values 
are given in Table VII . 

Figures 3.13 to 3.18 show the graphical results for 
carbon temperature, oxygen concentration and reaction rate 
for this case. 

TEMPERATURE SURFACE FROM GRAF3E 




Figure 3.13 Temperature vs X/L and time 
for thickness = 0.75 in 
Initial carbon temperature = 900 F. 
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TEMPERATURE SURFACE FROM GRHE3E 




Figure 3.14 Temperature vs X/L and time 
for thickness = 0.75 in 
Initial carbon temperature = 890 F. 

OXYGEN CONC. SURFACE FROM GRAF3E 




Figure 3.15 Oxygen concentration vs X/L and time 
for thickness = 0.75 in 
Initial carbon temperature = 900 F. 
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OXYGEN CONC. SURFRCE FROM GRRF3E 




Figure 3.16 Oxygen concentrat ion vs X/L and time 
for thickness = 0.75 in 
Initial carbon temperature = 890 F. 

RERCT10N RRTE SURFRCE FROM GRRF3E 




Figure 3.17 Reaction rate vs X/L and time 
for thickness = 0.75 in 
Initial carbon temperature = 900 F. 



62 



REACTION RATE SUREACE FROM GRAF3E 




Figure 3.18 Reaction rate vs X/L and time 
for thickness = 0.75 in 
Initial carbon temperature = 890 F. 

4 . CASE III-4 Thickness = 1.00 inches 

The program was run with several values of initial 
uniform graphite temperature . The temperatures bounding 
combustion and extinction are: 



for combustion Tc = 850 F 
for extinction Te = 840 F 

The variable parameters of this case attain maximum 
and minimum values as given in Table VIII . 

In this case, the graphical results for carbon 
temperature, oxygen concentration and reaction rate are 
shown in Figures 3.19 to 3.24 . 
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temperature surfrce frcf grp~3E 




Figure 3.19 Temperature vs X/L and time 
for thickness = 1.00 in 
Initial carbon temperature = 850 F. 



7E“ = ERR7URE SURFACE FRC" 1 GR pr 3E 




Figure 3.20 Temperature vs X/L and time 
for thickness = 1.00 in 
Initial carbon temperature = 840 F. 
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OXYGEN CONG. SURFACE PROA GRFF3C 




Figure 3.21 Oxygen concentration vs X/L and time 
for thickness = 1.00 in 
Initial carbon temperature = 850 F. 



CaYGEN conc. surface 




Figure 3.22 Oxygen concentration vs X/L and time 
for thickness = 1.00 in 
Initial carbon temperature = 840 F. 
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REACTION RATE SURFACE FROM GRAF3E 




3.23 Reaction rate vs X / L and time 
for thickness = 1.00 in 
Initial carbon temperature ’= 850 F. 

REACTION RATE SURFACE FROM GR a f"3E 




Figure 3.24 Reaction rate vs X/L and time 
for thickness = 1.00 in 
Initial carbon temperature = 840 F. 
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5. 



CASE III-5 Thickness = 2.00 inches 



The combustion and extinction process, in this case, 
are characterized by the following uniform initial graphite 
temperatures : 

for combustion Tc = 760 F 
for extinction Te = 750 F 

The graphite temperature and other parameters varied 
during these processes . The variation of these parameters 
characterize the system. The maximum and minimum values of 
these parameters are in Table IX . 

The graphical results for carbon temperature, oxygen 
concentration and reaction rate are shown in Figures 3.25 to 
3.30 . 

TEMPERATURE SURFACE FROM GRPF3E 



t 

5 *’• 

S'* 
I “*> 

I 

{ •« 




Figure 3.25 Temperature vs X/L and time 
for thickness = 2.00 in 
Initial carbon temperature = 760 F. 
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TEMPERATURE SURFACE FROM GRAF3E 




Figure 3.26 Temperature vs X/L and time 
for thickness = 2.00 in 
Initial carbon temperature = 750 F. 

OXYGEN CONC. SURFACE FROM GRAF3E 




Figure 3.27 Oxygen concentration vs X/L and time 
for thickness = 2.00 in 
Initial carbon temperature = 760 F. 
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OXYGEN CONC. SURFACE FROM GRAF3E 




Figure 3.28 Oxygen concentration vs X/L and time 
for thickness = 2.00 in 
Initial carbon temperature = 750 F. 

REACTION RATE SURFACE FROM GRAjcfc&E 




Figure 3.29 Reaction rate vs X/L and time 
for thickness = 2.00 in 
Initial carbon temperature = 760 F. 
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REACTION RATE SURFACE FROM GRAF3E 




Figure 3.30 Reaction rate vs X/L and time 
for thickness = 2.00 in 
Initial carbon temperature = 750 F. 

6 . CASE III-6 Thickness = 4.00 inches 

The minimum value of the initial graphite tempera- 
ture to start combustion of the porous medium and the 
maximum value of the initial graphite temperature that 
results in extinction for this case are: 



for combustion Tc = 690 F 
for extinction Te = 680 F 



Other variable parameters of this case are shown in 
Table X . The values of these parameters are characteristic 
of the combustion and extinction processes for this value of 
thickness. 

In this case, the graphical results show the 
behavior of carbon temperature , oxygen concentration and 
reaction rate in Figures 3.31 to 3.36 
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uRc SURFflCI FRO' - ' 



'PP 




Figure 3.31 Temperature vs 
for thickness = 4.00 
Initial carbon temperature 



X/L and time 
in 

= 690 F. 



:"°f==-^?p 



CP 



ir 




Figure 3.32 Temperature vs X/L and time 
for thickness = 4.00 in 
Initial carbon temperature = 680 F. 
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Figure 



Figure 



:rr:-Lr, 



SUPFrCE FROM GR C F3E 




j < 33 Oxygen concentration vs X/L and time 
for thickness = A . 00 in 
Initial carbon temperature = 690 F. 



G'T'EIb :Cr<C 










3.34 Oxygen concentration 
for thickness = 4.00 
Initial carbon temperature 



vs X/L and 
in 

= 680 F. 



time 
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PEP'CTIuN PATE 5 u =— CE FRO- 6RJ?)F3E 




Figure 3.35 Reaction rate vs X/L and time 
for thickness = 4.00 in 
Initial, carbon temperature .= 690 F. 

REACTION RATE SURFACE FROM GRAF3E 




Figure 3.36 Reaction rate vs 
for thickness = 4.00 
Initial carbon temperature 



X/L and time 
in 

= 680 F. 
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CASE III-7 Thickness 



6.00 inches 



7 . 

Below is given the initial uniform carbon tempera- 
ture for combustion and extinction: 

for combustion Tc = 650 seconds 
for extinction Te = 640 seconds 

The variable parameters for this case are given in 
Table XI 

The Figures 3.37 to 3.42 show the development of 
carbon temperature, oxygen concentration and reaction rate 
at each position during the combustion and extinction 
processes . 




Figure 3.37 Temperature vs X/L and time 
for thickness = 6.00 in 
Initial carbon temperature = 650 F. 
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TEMPERATURE SURPRCE PROM GRRF3E 




Figure 3.38 Temperature vs X/L and time 
for thickness = 6.00 in 
Initial carbon temperature = 640 F. 

cxy:iN sjpppcl r -c'i 



F igure 




3.39 Oxygen concentration vs X/L and time 
for thickness = 6.00 in 
Initial carbon temperature = 650 F . 
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OXYGEN CONC. SURFACE FROM GRAF3E 




Figure 3.40 Oxygen concentration vs X/L and time 
for thickness = 6.00 in 
Initial carbon temperature = 640 F. 





Figure 3.41 Reaction rate vs X/L and time 
for thickness = 6.00 in 
Initial carbon temperature = 650 F. 
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RERCT10N RRTE 3URFRCE FRCfl GPRF2E 




Figure 3.42 Reaction rate vs X/L and time 
for thickness = 6.00 in 
Initial carbon temperature = 640 F. 



D . SUMMARY 

Here, some observations are made about the effects of 
thickness on the combustion process. 

1 . Power Relation 

For each case a pair of values (thickness, combus- 
tion temperature) was obtained. Four plots of this data are 
shown in Figures 3.43 to 3.46 . 

Figures 3.43 and 3.44 show thickness versus combus- 
tion temperature , and Figures 3.45 and 3.46 show thickness 
versus the time it takes for the oxygen concentration at X/L 
= 1. to reach. 

It is observed that the relation between thickness 
and combustion temperature , figure 3.43 has a form which 
approximates an exponentially decreasing function. 

The log-log plot ( Figure 3.44 ) is close to a 

linear form, which yields the approximate power relation. 

The temperature is in degrees fahrenheit and thick- 
ness is in inches . 
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-0 16554 



(3.1) 



T = 850 L 

A comparison of the computer program’s results and 
the formula's results are given in Table III. The maximum 
difference between the analytical results and the power 
equation results is less than 3 percent over the range of 
thickness between 0.25 and 6.00 inches . 

This shows that the formula is a good approximation 
for these cases. 

2 . Combustion Start Time 

During the initial phase of the combustion process, 
the oxygen concentration at all points of the porous medium 
continuously decreases to zero, for all cases, the point of 
the porous medium that first reached zero was X/L = 1. The 

time it taken for that the oxygen concentration at X/L = 1. 
to reach zero is important because it is related to the time 
that combustion starts. 

Figure 3.45 , shows that this time does not increase 
linearly with the thickness. The log- log plot of these 
values, Figure 3.46 , shows that this log-log plot is 
approximately a straight line. 

The equation for this line is : 

T* = 226.8 L 1 . 683 (3.2) 

This equation can be used to give an estimate for 
the time it takes for combustion to start. 
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COMBUSTION TEMPERATURE (F) 

600 700 89O 000 1000 1100 1200 



THICKNESS VERSUS TEMPERATURE 




THICKNESS (INCHES) 



Figure 3.43 Thickness versus combustion temperature 

Retangular plot . 
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THICKNESS X TEMPERATURE 




Figure 3.44 Thickness versus combustion temperature 

Log- log plot . 
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TIME THAT OXYGEN REACH ZERO (S) 

1600 MOO *800 3000 MOO 



THICKNESS VERSUS T* 




Figure 3.45 Thickness vs time for oxygen concentration 

at X/L = 1. to reach zero 
Retangular plot. 
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TIME THAT OXYGEN REACH ZERO (S) 

i?o 1000 



THICKNESS VERSUS T* 




Figure 3.46 



Thickness vs 
reach zero 



time for oxygen 
log-log plot. 



concentration 
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IV. PERMEABILITY 



A. INTRODUCTION 

The permeability or hydraulic conductivity of a porous 
medium is the capacity of air flow through the porous 
medium. The air provides the oxygen for combustion . On the 
other hand , air flow provides for convective heat transfer 
within the medium . Thus air flow provides for heat genera- 
tion from combustion , as well as heat transfer by convec- 
tion . The interaction of these two mechanisms determines 
whether combustion or extinction will occur . If the heat 
generation of combustion is greater than the heat transfer 
of convection , then combustion will occur . If the heat 
transfer mechanism dominates , then extinction will occur . 
The equation used for permeability in the program is : 



Where m is the permeability , p is the porosity and z is the 
specific internal area . The expressions for porosity and 
specific internal area are : 



m= 0 . 2 p 3 / z 2 



(4.1) 



P = 1 - 7T / 6 (d/D ) 3 



(4.2) 



z = 1/2 7T (d/D 2 ) 



(4.3) 



z = 1/2 v (d 2 /D 3 ) 



(4.4) 



The equation 4.3 is for spherical particles and equation 4.4 
is for cylindrical particles . Equation 4.1 is the Kozeny 



equation for permeability of porous medium . This equation 
fails for highly porous fibrous media . 




Figure 4.1 Geometry of a typical cell. 

This equation fails for highly porous fibrous medias . 

The tortuosity ,r , is the ratio of the length of the 
flow path for a fluid particle to the straight line 
distance. The tortuosity depends on the ratio d/D .In this 
program tortuosity was assumed to have a constant value of 
1.4 . 

Normally the interior geometry of a porous medium can be 
inregular and complex . In such cases the internal structure 
is difficult to discribe analytically . The program used in 
this investigation models the porous medium as consisting of 
cylindrical fibers or spherical particles in a regular peri- 
odic geometry as shown in Figures 4.1 . 

During the combustion process carbon is being continu- 
ously consumed and as a result, the internal geometry is 
changing. The medium becomes more porous as the carbon 
decreases. Therefore, all properties which depend upon fiber 
or particle diameter are functions of time and position. The 
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model used in the program assumed that the carbon matrix 
remains rigid as the particles diameter decreases , and thus 
, the porosity increases with combustion . 

B . PROCEDURE 

To study the effects of permeability on combustion , all 
parameters of the program were maintained constant , except 
the permeability . In order to change the permeability and 
keep the porosity constant , it is necessary to keep the 
(d/D) ratio constant as seen in the equations 4.5 and 4.6 . 

This is accomplished by setting d/D equal toct • Then the 
numerator of equation 4.1 remains constant , however the 
denominator changes because d =a D and Z gives : 

Z = 1/2 Tr(a d/D 2 ) = constant / D (4.5) 

Z = 1/2 d 2 /D 2 )= constant / D (4.6) 

Thus , the dimensions of d and D were changed propor- 
tionally so that the initial ratio (d/D) remained constant 
in all cases . Thus the porosity remained constant while the 
permeability was different for each case . 

The cases studied in this section are : 



CASE IV- 1 
CASE IV- 2 
CASE IV- 3 
CASE IV-4 



PERMEABILITY 

PERMEABILITY 

PERMEABILITY 

PERMEABILITY 



0.00015 ft 2 (d=D=0 .0025 in) 
0.00058 ft 2 (d=D=0 . 0050 in) 
0.00233 ft 2 (d=D=0 .0100 in) 
0.00933 ft 2 (d=D=0 . 0200 in) 



For each values of permeability the program was run with 
different uniform carbon temperatures for all points (X/L) . 
This procedure was repeated until we obtaind the minimum 
temperature that results in combustion . This temperature 
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is called the combustion temperature (T c ) 
occurs at temperatures below T c . 



Extinction 




(4.7) 



In this section this difference (5 ) was chosen equal at 
10 F then : 



C. RESULTS 

For each case , the program provided the values of the 
carbon and air temperatures inside the porous medium , as 
well as other parameters which varied during the combustion 
process . 

1 . CASE IV-1 Permeability = 0.00015 ft^ (d=D=0.0025 in) 
For this value of permeability , the maximum carbon 
temperature that does not result in combustion (extinction 
temperature) and the minimum temperature that results in 
combustion (combustion temperature) are : 



4.7 . These graphical surfaces show the carbon temperature , 
oxygen concentration and reaction rate at each position and 
time . The form of these surfaces shows the different 
behavior of these three parameters during the extinction and 
combustion processes. 



T c = T e + 10 (F) 



(4.8) 



for combustion T c = 760 F 

for extinction T 0 = 750 F 



The graphical results are shown in Figures 4.2 to 
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TEMPERATURE SURFACE FROM GRAF3E 




Figure 4.2 Temperature vs X/L and time 
for permeability = 0.00015 ft^ 
initial carbon temperatiure = 760 F. 



TEMPERATURE SURFACE FROM GRAF3E 




Figure 4.3 Temperature vs X/L and time 
for permeability = 0.00015 ft^ 
initial carbon temperatiure = 750 F. 
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OXYGEN CONC. SliRFRCE FROM GRRF3E 




Figure 4.4 Oxygen concentration vs X/L and time 
for permeability = 0.00015 ft“ 
initial carbon temperatiure = 760 F. 

OXYGEN CONC. SURFRCE FROM GRPF3E 




Figure 4.5 Oxygen comcentrat ion vs X/L and time 
for permeability = 0.00015 ft^ 
initial carbon temperatiure = 750 F. 
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REACTION RATE SURFACE FROM GRPF3E 




Figure 4.6 Reaction rate vs X/L and time 
for permeability = 0.00015 ft 
initial carbon temperatiure = 760 F. 

REACTION RATE SURFACE FROM GRAF3E 




Figure 4.7 Reaction rate vs X/L and time 
for permeability = 0.00015 ft 
initial carbon temperatiure = 750 F. 
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2. CASE IV- 2 Permeability = 0.00058 ft 2 (d=D=0.0050 in) 

O 

Here the value of permeability is 0.00058 ft which 
results from a pore diameter of 0.0050 in . For this case 
the combustion and extinction temperatures are : 



for combustion T c - 850 0 F 

for extinction T 0 = 840 0 F 

From the numerical output of the program , three 
graphical surface s were created ; the carbon temperature , 
oxygen concent rat ion and reaction rate as functions of time 
and position . These surfaces provide information of the 
development of these three dependent variables during the 
combustion and extinction process . These surfaces are shown 
on Figures 4.8 to 4.13 . 

TEIMPERRTURE SURFRCE from GRRF3E 




Figure 4.8 Temperature vs X/L and time 
for permeability = 0.00058 ft^ 
initial carbon temperatiure = 850 F. 
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TEMPERATURE SURFACE FROM GRA r 3E 




Figure 4.9 Temperature vs X/L and time 
for permeability = 0.00058 ft 
initial carbon temperatiure = 840 F. 

OXYGEN CONC. SURFACE FROM GRA r 3E 




Figure 4.10 Oxygen concentration vs X/L and time 
for permeability = 0.00058 ft^ 
initial carbon temperatiure = 850 F. 
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OXYGEN CONC. SURFACE FROM GRAF3E 




Figure 4.11 Oxygen comcentration vs X/L and time 
for permeability = 0.00058 ft^ 
initial carbon temperatiure = 840 F. 



REACTION RATE SURFACE "ROM GRA r 3E 




Figure 4.12 Reaction rate vs X/L and time 
for permeability = 0.00058 ft^ 
initial carbon temperatiure = 850 F. 
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REACTION RATE SURFACE FROM GRAF3E 




Figure 4.13 Reaction rate vs X/L and time 
for permeability = 0.00058 ft^ 
initial carbon temperatiure = 840 F. 

3 . CASE IV-3 Permeability = 0.00233 ft 2 (d=D=0.0100 in) 
For this case the initial permeability is 0.00233 

O 

ft . With this value of permeability the combustion and 
extinction temperatures are : 



for combustion T c = 980 0.F 
for extinction T £ = 970 0.F 

The variation of the dependent variables during the 
combustion and extinction cases is shown by a numerical 
output of the program . These results are also shown in 
three graphical surfaces ; carbon temperature , oxygen 
concentration and reaction rate on each position (X/L) and 
time . These graphical results are shown in Figures 4.14 to 
4.19 
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TEMPERATURE SURFACE FROM 3RAF3E 




Figure 4.14 Temperature vs X/L and time 
for permeability = 0.00233 ft^ 
initial carbon temperatiure = 980 F. 

TEMPERATURE SURFACE FRC^ GRAF3E 




Figure 4.15 Temperature vs X/L and time 
for permeability = 0.00233 ft^ 
initial carbon temperatiure = 970 F. 
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QX r GEN CONC. SURFRCE ^RC^ GR Q F3E 




Figure 4.16 Oxygen concentration vs X/L and time 
for permeability = 0.00233 ft^ 
initial carbon temperatiure = 980 F . 

OXYGEN CQNC . 5URFPCE r RCT GR C F3E 




Figure 4.17 Oxygen comcentration vs X/L and time 
for permeability = 0.00233 ft^ 
initial carbon temperatiure = 970 F. 
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PERCH ON RATE SUR ro CE FRCM G^AFSE 




Figure 4.18 Reaction rate vs X/L and time 
for permeability = 0.00233 ft 
initial carbon temperatiure = 980 F. 

REACTION RATE SURFACE FRCF GRA r 3E 




Figure 4.19 Reaction rate vs X/L and time 
for permeability = 0.00233 ft 
initial carbon temperatiure = 970 F. 
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4. CASE IV -4 Permeability = 0.00933 ft 2 (d=D=0.0200 in) 

The combustion and extinction temperatures for this 

2 

case , with initial permeability equal to 0.00933 ft , are: 

for combustion T c = 1180 0.F 

for extinction T e = 1170 0.F 

for this case , it was not possible to obtain the graphical 
surfaces due to limitations of the graphical program . 

D . SUMMARY 

Here , some observations about the effects of perme- 
ability on combustion can be understood by looking at all 
the results at one time . Table XIII contains the results of 
all cases . 

The Table lists the temperatures for extinction and 
combustion , the permeability (m) , the dimension of the unit 
cell (d=D) and time t* . The t* time is the time at which 
the oxygen concentration at a point reaches zero first . 
That point of the porous medium that first reaches zero 
oxygen concentration in all these cases is the the position 
point X/L = 1 . 

1 . Power Relation 

The pair of values (permeability , combustion 
temperature) for each case was plotted on a cartesian graph. 
An observation of this curve shows that the combustion 
temperature increases with increasing permeability. Figure 
4.20 shows that the combustion temperature increases mono- 
tonically and approaches an assindoticptate as permeability 
increases . The plot of these same points on a log-log 
graph , Figures 4.21 , shows that its form is not quite a 

straight line . 

It is observed that the relation between combustion 
temperature (T c ) and permeability (m) yields the approximate 
power relation (equation 4.9) : 
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(4.9) 



T c = 1880 m 0 . 105 

A comparation of the values of the combustion 
temperature obtained by the program and by the power rela- 
tion is presented in Table XIV . This Table shows that the 
difference between these results fail between 1.2 to 3.2 
percent . This result shows that the above power relation is 
a good approximation to estimate the temperature which will 
start combustion for a given permeability . 

2 . Combustion Speed 

In this section , the question of how permeability 
affects the time it takes for combustion to begin is 
addressed . We meassure the rapidilly of the combustion 
process by the t* time previously defined . Recall that t* 
is the time at which the oxygen first reaches zero at some 
point in the porous medium . In all cases the position X/L = 
1 is the first point to achieve zero oxygen concentration . 
The plots of permeability (in) versus t* , Figures 4.22 and 
4.23 , show that t* decays exponentially with increasing 
permeability and approaches an assymptotic value of t* = 0 
as m goes to infinity . 

A very interesting observation is obtained from the 
results of Table XIII . It is noted that the m^+^/m^ (for i 
= 1,2,3) remains pratically constant (3.84 , 4.01 , 4.00) 

while the t*^ + ^/t*^ ratio takes on the variable values 
1/ (i+1) (1/1.91 , 1/2.65 , 1/3.56) . 

An attempt was made to obtain an additional run with 

O 

the permeability equal to 0.03733 ft , however the combus- 
tion activity was too rapid to achieve meaningful results . 
Combustion was achieved so rapidly that the gradient were 
too large for the program to handle . 

The permeability is a parameter that has great 
influence on the combustion process for porous medium . 
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PERMEABILITY VERSUS TEMPERATURE 




Figure 4.20 Permeability versus combustion temperature 

retangular plot. 
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COMBUSTION TEMPERATURE (F) 



PERMEABILITY VERSUS TEMPERATURE 




Figure 4.21 Permeability versus combustion temperature 

log log plot . 
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TIME THAT OXYGEN REACH ZERO (S) 

IM fOO 800 800 800 400 400 



PERMEABILITY VERSUS T* 




0.000 0.008 0.004 0.000 0.000 0.010 



PERMEABILITY (FT.SQ) 

Figure 4.22 Permeability versus t* 
Retangular plot. 
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PERMEABILITY VERSUS T* 



o 




Figure 4.23 Permeability versus t - 
Log-log plot. 
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V. REACTION ORDER 



A. INTRODUCTION 

The expression for reaction rate for carbon reacting in 
air used in the program was the expression resulting from 
the experimental data of Parker and Hotel . 



Reaction order is the exponent of the partial pressure 
of oxygen in the expression of the reaction rate . The 
program permits the user to select any value for this param- 
eter . In all of the previous chapters , a value of n equal 
to 0.5 was used . In this chapter we investigate the effect 
of n on the combustion process . 

B. PROCEDURE 

In this section , all parameters are fixed to their 
basic values and only the reaction order is changed in each 
case . 

The program was run , for each value of reaction order , 
with different value of the uniform initial carbon tempera- 
ture until we obtained the minimum temperature which results 
in combustion . This temperature is the combustion tempera- 
ture (T c ) . The maximum temperature that results in extinc- 
tion is the extinction temperature (T e ). they differ by an 
infinitesimals (T c = T e +5) • The magnitude of<5 was chosen 
equal to 10 F for this study . 

In order to study the influence of the reaction order on 
the combustion process , four different values of this 
parameter were chosen . These cases are : 




(5.1) 
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CASE V-l 
CASE V- 2 
CASE V- 3 
CASE V-4 



Reaction Order = 0.25 
Reaction Order = 0.50 
Reaction Order = 0.75 
Reaction Order = 1.00 



The fixed parameters used in all cases are : 

Ambient temperature =80 F 
Ambient pressure = 2117 lb/ft. sq 
Tortuosty = 1.400 
Filamento diameter = 0.0004167 ft 
Thickness of matrix laminette = 0.0004167 ft 
Thickness of porous medium = 0.02083 ft 
Gas constant for air = 53.34 lbf ft/lbm R 
Conductivity of filament = 86 BTU/lbm H.F. 

Specific heat of filament = 0.703 lbm/Cf 
Emissivity of filament = 0.90 
Shape factor for internal HFxFER =1.00 
Characteristc length of medium = 1 inchess 
Heat of reaction = 14090 BTU/Cf 
Stochiometric ratio (fuel/air) = 0.375 
Reaction coefficient = 2065000 lbm/Cf H 
Activetion energy coefficient = 28840 deg R 
Pressure differential across thickness = 50 lb/ft. sq 
Initial uniform oxygen concentration = 0.172 lbm/ft 

C. RESULTS 

For each cases the results were obtained in both numer- 
ical and graphical form . These results show the variation 
of some parameters during the combustion and extinction 
processes . The results obtained for each value of reaction 
order are presented below 

]. . CASE V-l Reaction Order = 0.25 

With this value of reaction order , the combustion 
and extinction processes define the following characteristic 
temperatures : 
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980 o.F 
970 o.F 



for combustion T c = 
for extinction T e = 

The analysis of the numerical output of the program 
shows that some parameters varied during the combustion and 
extinction processes . The limits of these parameters show 
the character of the particular kind of process which 
occured . 

The graphical results show three surfaces that 
represent the evolution of the carbon temperature , oxygen 
concentret ion and reaction rate at each position during the 
transient history . These graphical surfaces are shown in 
Figures 5.1 to 5.6 . 

TEMPERRTURE 5UR r PCE FRQr 3R cr 3E 




Figure 5.1 Temperature vs X/L and time 
for reaction order = 0.25 
Initial carbon temperature = 980 F. 
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TEMPERATURE SURFACE FROM GRAF3E 




Figure 5.2 Temperature vs X/L and time 
for reaction order = 0.25 
Initial carbon temperature = 970 F. 

OX T GEN CONC. SURFACE ~RC M . C-RRF3E 




Figure 5.3 Oxygen concentration vs X/L and time 
for reaction order = 0.25 
Initial carbon temperature = 980 F. 
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OXYGEN CONC. SURFACE FROM GRAF3E 




Figure 5.4 Oxygen concentration vs X/L and time 
for reaction order = 0.25 
Initial carbon temperature = 970 F. 



REACTION PATE SURFACE ARC” GRAF3E 




Figure 5.5 Reaction rate vs X/L and time 
for reaction order = 0.25 
Initial carbon temperature = 980 F. 
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REACTION RATE SURFACE FROM GRAF3E 




Figuie 5.6 Reaction rate vs X/L and time 
for reaction order = 0.25 
Initial carbon temperature = 970 F. 



2 . CASE V-2 Reaction Orde r = 0.50 

For reaction order equal to 0.50 , the initial 
uniform carbon temperatures that separate the combustion and 
extinction processes are : 



for combustion T c = 980 F 
for extinction T e = 970 F 

The program’s numerical output included the values 
of some parameters that varied during the transient process. 
The variation of the values characterize the combustion or 
extinction process. The graphical results are shown in 
three surfaces: carbon temperature, oxygen concentration and 
reaction rate versus position (X/L) and time. These graph- 
ical surfaces are shown in figures 5.7 to 5.12 . 
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TEMPERATURE SURFACE FROM GRAF3E 




Figure 5.7 Temperature vs X/L and time 
for reaction order = 0.50 
Initial carbon temperature = 850 F. 

TEMPERATURE SURFACE FROM GRAF3E 




Figure 5.8 Temperature vs X/L and time 
for reaction order = 0.50 
Initial carbon temperature = 840 F. 
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OXYGEN CONC. SURFACE FROM GRAF3E 




Figure 5.9 Oxygen concentration vs X/L and time 
for reaction order = 0.50 
Initial carbon temperature = 850 F. 

O/rGEN CONC. SURFACE FROM GRRF3E 




Figure 5.10 Oxygen concentration vs X/L and time 
for reaction order = 0.50 
Initial carbon temperature = 840 F. 
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REACTION RATE SURFACE FROM GRAF3E 




Figure 5.11 Reaction rate vs X/L and time 
for reaction order = 0.50 
Initial carbon temperature = 850 F. 

REACTION RATE SURFACE FROM GRAF3E 




Figure 5.12 Reaction rate vs X/L and time 
for reaction order = 0.50 
Initial carbon temperature = 840 F. 
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3. CASE V-3 Reaction Order = 0.75 



The reaction order was set equal to 0.75 for this 
case. The combustion and extinction temperatures are: 



for combustion Tc = 750 F 
for extinction Te = 740 F 



A visualization of the development of the carbon 
temperature, oxygen concentration and reaction rate at each 
position and time are shown in the graphical surfaces of 
Figures 5.13 to 5.18 . 



TEMPERATURE SURFACE FROM GRAF3E 




Figure 5.13 Temperature vs X/L and time 
for reaction order = 0.75 
Initial carbon temperature = 750 F. 
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TEMPERRTURE 5URFPCE PROP GRRF3E 




Figure 5.14 Temperature vs X/L and time 
for reaction order = 0.75 
Initial carbon temperature = 740 F. 

OXYGEN CONC. SURFRCE FRQR GRRF3E 




Figure 5.15 Oxygen concentration vs X/L and time 
for reaction order = 0.75 
Initial carbon temperature = 750 F. 
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OXYGEN CONC. SURFACE from GRAF3E 




Figure 5.16 Oxygen concentration vs X/L and time 
for reaction order = 0.75 
Initial carbon temperature = 740 F. 

REACTION RATE SURFACE FROM GRAF3E 




Figure 5.17 Reaction rate vs X/L and time 
for reaction order = 0.75 
Initial carbon temperature = 750 F. 
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REACTION RATE SURFACE FROfl GRAF3E 




Figure 5.18 Reaction rate vs X/L and time 
for reaction order = 0.75 
Initial carbon temperature = 740 F. 

4. CASE V-4 Reaction Order = 1.00 

For this case, the temperature of the combustion and 
extinction processes are: 



for combustion T c = 670 F 
for extinction T £ = 660 F 

The Figures 5.19 to 5.24 show the behavior of the 
carbon temperature, oxygen concentration and reaction at 
each position and time. 

D . SUMMARY 

1 . Power Relation 

For each case, a pair of values (reaction order, 
combustion temperature) was obtained. The points were plotted 
in cartesian coordinates in Figure 5.25 . A log-log plot is 
shown in Figure 5.26 , and a semi log plot is shown in 

Figure 5.27 . It is observed that the relation between 
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TEMPERATURE SUREACE PROM GRAE3E 




Figure 5.19 Temperature vs X/L and time 
for reaction order = 1.00 
Initial carbon temperature = 670 F. 

reaction order (n) and combustion temperature (T c ) of figure 
5.27 yields the approximate power relation (equation 5.2) . 

Table 5.25 shows a resume of each case, combustion 
temperature, extinction temperature and the value of the 
combustion temperature obtained from Equation 5.2 and the 
percentage error . An observation of these results leads to 
the conclusion that this power relation is a good estimation 
for combustion temperature. 
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TEMPERATURE SURFACE FROM GRAF3E 




Figure 5.20 Temperature vs X/L and time 
for reaction order = 1.00 
Initial carbon temperature = 660 F. 

0XTGEN CONC. SURFACE FROM GRAF3E 




Figure 5.21 Oxygen concentrat ion vs X/L and time 
for reaction order = 1.00 
Initial carbon temperature = 670 F. 
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OXYGEN CONC. SURFACE FRCA GRA r 3E 




Figure 5.22 Oxygen concentration vs X/L and time 
for reaction order = 1.00 
Initial carbon temperature = 660 F. 

| 

REACTION RATE SURFACE FROfl GRAF3E 




Figure 5.23 Reaction rate vs X/L and time 
for reaction rate= 1.00 
Initial carbon temperature = 670 F. 
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REACTION RATE SURFACE FROM GRAF3E 




Figure 5.24 Reaction rate vs X/L and time 
for reaction order = 1.00 
Initial carbon temperature = 660 F. 



T c = 10 . 05-0 . 22n) ( p ) 



(5.2) 
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TABLE XV 



COMBUSTION AND EXTINCTION TEMPERATURE 
FOR REACTION ORDER CASES 

temperatures ( F) 





extinction 


combustion 


% error 


reaction 


program 


program 


power 




order 


output 


output 


relation 




0.25 


970 


980 


980 


0.0 


0.50 


840 


850 


863 


1.6 


0.75 


740 


750 


760 


1.4 


1.00 


670 


670 


670 


0.0 
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COMBUSTION TEMPERATURE (F) 

•00 MO 700 700 800 800 800 800 1000 



REACTION ORDER VERSUS TEMPERATURE 




I ^ """ — ■ : 

0.00 O-tS 0.00 0.70 1.00 IM 



REACTION ORDER 

Figure 5.25 Reaction order versus temperature 

Retangular plot. 
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COMBUSTION TEMPERATURE (F) 



REACTION ORDER VERSUS TEMPERATURE 




Figure 5.26 Reaction order versus temperature 

Log-log plot. 
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REACTION ORDER VERSUS TEMPERATURE 




Figure 5.27 Reaction order versus temperature 

Semi log plot. 
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VI. CONCLUSIONS 



In each of the previous chapters, the effects of a 
single parameter on the combustion behavior of a porous 
medium was investigated. In Chapter 2, we considered the 
effect of heat flux ( SQ ) and its duration (TQ) on 
combustion. In Chapter 3, the effect of medium thickness (L) 
on behavior was investigated. Chapter 4 was a study of the 
effect of permeability (in) on system behavior; and finally 
in Chapter 5, the effect of reaction order (n) was 

investigated . 

The major conclusions from these investigations are: 

- The relation between SQ and TQ is not linear , a power 
relation gives rasonable approximation. 

- During the combustion process, for the SQ cases, 

initially there is a large temperature gradient across 
the medium with T ( X/L = 0.) >> T ( X/L =1.) . With 
passing time, this gradient decreases until there is a 
uniform combustion temperature over the entire medium. 
When this occurs , there is no oxygen within the 

medium, and the medium is continuosly consumed by 
combustion. 

- The relation between thickness and combustion tempera 

ture is not linear . Increasing the thickness results 
in decreasing of combustion temperature. If the 

thickness increased with a constant rate c ( L^ + ^/l.£ = 
c ) the time that marks the start of combustion is not 
constant , it increases with increasing thickness ( 
i.e. t'i + l/t'i = b(L) ) where b depends upon the 

thickness . 

- The relation between permeability and minimum combus 
tion temperature is not linear. 

- An increase of the permeability shows that the combus 
tion temperature increases to an asymptotic limit. 
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- Permeability has a significant influence on the speed 
of the combustion process. When the permeability is 
increased, the time that marks the start of the 
combustion process decreases. For large values of 
permeability, the • combustion may occur in a fraction 
of a second. 

- If the permeability is increased with a constant rate 

c( = c m^ ) , the time that measures the start of 

combustion decreases . 

- The relation between reaction order (n) and combustion 
temperature (Tc) is not linear. 

- An observation of the graphical surfaces, for all 

cases of this investigation , shows a similar form for 
all combustion and extinction cases. For all 
combustion cases, it was noted that the oxygen within 
the porous medium being consumed with passing time. 
The first point at which the oxygen is depleted is 

always X/L = 1. With passing time the point with 

smaller X/L reach zero oxygen, until all points have 
zero oxygen. At this time the only oxygen is at the 

interface of the porous medium and the ambient 
environment. At this time the combustion process 
becomes a surface recession phenomena. At the same 
time the oxygen is being depleted the temperature is 
increasing. Eventually when all the interior oxygen 
has been consumed, there is a temperature gradient 
across the medium with T c (X/L=0) > T C (X/L=1). 

This thesis concludes with two recomendations . First , 
it is recommended that additional parameter studies he 

conducted . For example , it would be useful to determine 
the effects of pressure differential , conductivity , and 
other thermophysical properties on combustion . Another 

recommendation is to attempt to obtain a relation between 
combustion temperature and several parameters . 
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